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Abstract 

The radiative dileptonic decays B s {Bd) — > (I = e, //) are investigated within 

the Standard Model . The transition formfactors are calculated in the framework of 
the light cone QCD sum rules method and it is found that the branching ratios are 
B(B S -» e+e'j) = 2.35 x 10~ 9 , B(B S -> /i+/x _ 7) = 1-9 x 10 -9 , B(B d -» e+e~7) = 
1.5 x 10 -10 and B(Bd — ► /i + /i~7) = 1.2 x 10~ 10 . A comparison of our results with 
the constituent quark model predictions on the branching ratios is presented. 
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1 Introduction 



The Flavour Changing Neutral Current (FCNC) processes are the most promising field 
for testing the Standard Model (SM) predictions at loop level and for establishing new 
physics beyond that (for a review see fl[] and references therein) . At the same time the rare 
decays provide a direct and reliable tool for extracting information about the fundamental 
parameters of the Standard Model (SM), such as the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix elements V t d, V ts and V u b [0. 

After the experimental observation of the b — > 57 || and B — > X s 7 [|J processes, the 
interest is focused on the other possible rare 5-meson decays, which are expected to be 
observed at future £>-meson factories and fixed target machines. In addition to being used 
in the determination of the CKM matrix elements, the rare 5-meson decays could play 
an important role in extracting information about some hadronic parameters, such as the 
leptonic decay constants fs a and fs d . Pure leptonic decays B s — > ji + ^ and B s — > e + e~~ are 
not useful for this purpose, since these decays are helicity suppressed and as a result they 
have branching ratios B(B S — > /U + /i~) ~ 1.8 x 1CT 9 and B(B S — > e + e~) ~ 4.2 x 1CT 14 ||. 
For Bd meson case the situation becomes worse due to the smaller CKM angle. Although 
the process B s — > t + t~, whose branching ratio in the SM is B(B S — > t + t~) = 8 x 10~ 7 
H, is free of helicity suppression, its observability is expected to be compatible with the 
observability of the B s —>■ fi + n~ decay only when its efficiency is better than 10 -2 . 

When a photon is emitted in addition to the lepton pair, no helicity suppression exists 
anymore and larger branching ratios are expected. For that reason, the investigation of 
the B s u) — > / + /~7 decay becomes interesting. The branching ratios of these processes 
depend quadratically on the leptonic decay constants of B mesons and hence it could be a 
possible alternate in determining f Ba and fs d - In 0, these decays are investigated in the 
SM using the constituent quark approach and it is shown that the diagrams with a photon 
radiation from the light quark give the dominant contribution to the decay amplitude which 
is inversely proportional to the constituent light quark mass. However the concept of the 
"constituent quark mass" is itself poorly understood. Therefore, any prediction on the 
branching ratios, in the framework of the above mentioned approach, is strongly model 
dependent. 

In this work, we investigate the B s u\ — > / + / _ 7 processes practically in a model indepen- 
dent way, namely, within the framework of the light cone QCD sum rules method (more 
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about the method and its applications can be found in a recent review ||). The paper is 
organized as follows: In sect. 2 we give the relevant effective Hamiltonian for the b — ► ql + l~ 
decay. In sect. 3 we derive the sum rules for the transition formfactors. Sect. 4 is devoted to 
the numerical analysis of the formfactors, and the calculation of the differential and total 
widths for the B q — > (q = s, d) decays. In this section we also present a comparison 

of our results with those of |7j . 

2 Effective Hamiltonian 

The most important contribution to B q — > Z + Z~7 (I = e, //) stems from the effective Hamilto- 
nian which induces the pure leptonic process B q . The short distance contributions 
to b — > l + l~q decay, comes from the box, Z-boson and photon mediated diagrams (Fig.l). 



The QCD corrected quark level amplitude in the SM can be written as R [Hi 



V2n 



M = -^V tb V* 



Cl f f{q llx P L b)T llx l + C 10 q^P L bl lfMl5 l - 

(1) 



c 7 _ 

-2—qia^p v {m h P R + m q P L )bl'y li l . 



Here Pl(r) = [1 ~ (+)7s] /2 , and p is the momentum of the lepton pair. The analytic 
expressions for all Wilson coefficients can be found in || [TOj . In further considerations we 
shall neglect the mass of the light quarks. 

As we have already noted, the pure leptonic processes B q —>■ (I = e, fx) are helicity 
suppressed. If a photon is attached to any of the charged lines in Fig.l, the situation 
becomes different; helicity suppression is overcome. If a photon is emitted from the final 
charged lepton lines, it follows from the helicity arguments that the amplitude of such 
diagrams must be proportional to the lepton mass mi (I = e, /i). Therefore the contribution 
of such diagrams are negligible. When a photon is attached to any charged internal line, 
the contributions of these diagrams will be strongly suppressed by a factor of ml/m^y in the 
Wilson coefficients, since the resulting operators have dimension 8, which are two orders 
higher than usual operators in (1). So, we conclude that the main contribution comes 
from the diagrams in Fig.l with a photon radiation from the initial quark lines. Thus the 
corresponding matrix element for the process B s ^) / + / _ 7 can be written as 
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(l\M\B) 




(2) 



These transition amplitudes can be written in terms of two independent, gauge invariant 
(with respect to the electromagnetic field) structures: 



Here, e M and stand for the polarization vector and the momentum of the photon, p is 
the momentum of the lepton pair, g{p 2 ) , gi{p 2 ) , and f{p 2 ) , fi(p 2 ) describe the parity 
conserving and parity violating formfactors. Thus, the main problem is to calculate the 
formfactors g , g\ and / , f\ including their momentum dependence. For this aim we will 
employ the light cone QCD sum rules method. 

Note that the formfactors g and / are calculated in the light cone QCD sum rules in 
fill • Therefore we concentrate ourselves to the calculation of formfactors g\ and fi induced 
by the magneto-dipole interaction. 



3 Sum rules for the transition formfactors f\(p 2 ) and 



According to the QCD sum rules ideology, in order to calculate the transition formfactors 
fi{p 2 ) and gi(p 2 ), it is necessary to write the representation of a suitable correlator function 
in the hadronic and quark-gluon languages. We start by considering the following correlator 
function 




(3) 
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This correlator can be calculated in two different ways. On one side we insert to H^(p, q) 
the hadronic states with B meson quantum numbers. Then we have 

m B Ib 1 

n M (p,g) = ^ m ^ _^ + ^ 2 ('y(q)\qiir lia p a (l + ^)b\B{p + q)) 

c m B q fB q 1 x 

m b rn 2 Bq - (p + q) 2 
x < e^ a e a pf3q a — + +i e^pq) - (e p)q^ — \ . (5) 



m B 

In deriving eq(5) we used 

TB a 



m B 



™bJi 



(B\bt l5 q\0) = 

m b 

On the other hand, the correlation function (4), can be calculated in QCD at large Eu- 
clidean momenta (p + q) 2 . In general, the correlator (4) can be decomposed into the parity 
conserving and parity violating parts 

n M (p, q) = e mf3 xe* a ppq>jli + i [e* (pq) - q^(e*p)] IT 2 (6) 

Equating eqs.(5) and (6) we get sum rules for the formfactors gi(p 2 ) and fi(p 2 )- 

Let us start calculating U^(p,q) from QCD side. The virtuality of the heavy quark 
in the correlator function under consideration, is large and of order m\ — (p + q) 2 . Thus, 
one can use the perturbative expansion of the heavy quark propagator in the external field 
of slowly varying fluctuations inside the photon. The leading contribution is obtained by 
using the free heavy quark propagator in eq.(4). Then we have 

H-n(p,q) = J ( 2?r )4 ( m 2 _ fc2 ) (7lg(a)KWa(l + TsXff + m b )i-f 5 q(0) |0) 

rd i xd A k e l ( p ~ k)x \ 
= -J (^pyPal^l^^l + 75)?(0)|0> - 

-k p (l\q{x)Wp(l- 75)9(0) |0) I (7) 

In this equation a path ordered gauge factor between the quark fields is omitted, since in 
the Fock-Schwinger gauge x^A^(x) = 0, where A^(x) is the external electromagnetic field, 
it is irrelevant. 

The diagrams (a) and (b) in Fig. 2 describe only the short distance (perturbative) part 
of these matrix elements corresponding to the photon emission from the freely propagating 
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heavy and light quarks. The non-perturbative contributions correspond to the propagation 
of the light quark in the presence of external electromagnetic field (Fig. 2c and 2d). 

We consider now the perturbative contributions. For the diagrams (2-a and 2-b) we can 
write down the double dispersion representation 

ds dt p\ ' '(s, t) 



n (l,2) 



+ subtr. terms. 



(8) 



s — (p + q) 2 } (t - p 2 

Here, superscripts 1, and 2 correspond to the contributions of the spectral densities to the 
structures e IMJl px e ^Pi3Qx and e* (pq) — q^(e*p) respectively. 



For calculating the spectral densities p\ and pn we use the method given in [pL2| . After 
a rigorous calculation for spectral densities, we have 



16n 



-ee q s5{t — s) 



(9) 



2Nc X(4 \(, ^ 

-^eebSdit - s) 1 



(10) 



P?\s,t) 



2iVc L(+ \d m < 
—ee q l5(t-s) 1-- 



-8'{t 



mi 



(s - ml) s 



(11) 



P { n\s,t) 



2N C 



eeA 5{t — s) 



mi 



3 1 

—s H — mi ) — 2m? In 
2 2 b ' b 



mi 



-5'(t 



mi 



(s 2 + smfj 



2sm b In 



mi 



(12) 



In eqs.(9-12) pi and pu corresponds to the interaction of the photon with the light and b 
quarks, iV c = 3 is the color factor, e q and the electric charge of the light and b quarks 
and is the mass of the 6-quark 5'{t — s) = 4:3 (t — s). 

Next consider the non-perturbative contributions. From eq.(7) it follows that the non- 
perturbative contributions are expressed via the matrix elements of the gauge invariant 
nonlocal operators, sandwiched in between the vacuum and the photon state. These matrix 



elements define the following light cone photon wave functions ( [jlQ| , P^3| , see also the first 
reference in ||11|| ): 
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(l\q(xW a q(0)\0) = iee q {qq) f due iu(lx 

J 

{(e M g a - e a q^) [x4>{u] + x 2 (g 1 (u) - g 2 (u))\ + 

+g 2 (u) [qx(e^x a - e a x^) + ex(x^q a - x a q^)} } and 
1 r 1 

(l\q{xh^ 5 q{0) |0) = -ee^ x e aPp x x f j g±(u)e mqx . (13) 

Here x is magnetic susceptibility of the quark condensate, 4>{u) , g±(u) are the leading twist 
t = 2 photon wave functions, g±(u) and #2 (w) are the two particle r = 4 wave functions. 
Note that for calculating the matrix elements 

(7(g)|g(x)cr MQ 7 p (l - 7s)g(0)|0) and 
(7(g)l?(aO<^75<z(0)|0) 

we use the following identities: 

<V*75 = -itfiaXpVXp , (14) 

<V*7p = i{l»9ap - lag^p) + e MapA 7A75 • (15) 

After lengthy calculations for n x and il 2 we get the following results, which describe 
the non-perturbative contributions: 



III = m b ee q {qq) [ du\ - + 8m 2 



9i(u)-g 2 (u) 4(m 2 b -p 2 ^ 



A 3 



A 3 



-92 



1 r + mj 

A A 2 



no = 



r Jo 

m b ee q (qq) dul - X ^ 



9±(u) + ee b m b (qq) 



(m 2 b - p 2 ) [m 2 b - (p + <?) 2 ] 



, (16) 



+ 8m; 



4(m 2 +p 2 ) 



A 3 



A 3 



"#2 



2 2-upg 
A 



+ ee b m b (qq)- 



(17) 



(m 2 -p 2 ) [m 2 - {p + q) 2 ] 

Here A = m 2 — (p + -ug) 2 . Last term in eqs.(16) and (17) describes the case when a 
photon is emitted from the heavy quark (see Fig. 2d). Collecting eqs.(9-12) and (16-17) we 
finally get the following expressions for the invariant functions Hi and n 2 : 
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Ill 



N c e r du ™l-fu / 



16vr 2 



mtu 



u 2 A\ 



ml — p 2 u j 



e 1 + 



mtu 



m\ — p 2 u j 



lei 



+ m b ee q (qq) / — 



du f l0 25 , i(«) -P2(m) (m 2 b -p 2 ) 



e r du 
V J ~u 



1 p 2 + ml 
~K~x + wA 2 



Ai +8m& « 2 A? 

+ ee b m b (qq)- 



-4- 



u 2 A? 
1 



-5-2 



m 2 - p 2 ) [m 2 - (p + g) 2 ] 



18) 



n. 



N c e ri 



167T 2 7o Ai (m 2 — p 2 ) 



' m 2 b -p 2 u 2 N 



■umr 



m 2 — p 2 u / 



-,2„-, 



+ e a 



mi — p u l mi — p u 



u 



+ Ae h mlp 2 ln ( m& J' - 
V m 6 M , 



+ 



+ m b ee q (qq) / — , 

JO W I Ai 



cfuf , _ 2 gx(u)-g 2 (u) (m 2 b +p 2 ) 



+ 8m 



2 A? 



2 A? 



#2 



3p 2 + 



e r 1 du 
V Jo ~ 



2 | 2(pq) 
Ai mA? 



^jl(u) + ee b m b (qq) 



(ml - p 2 ) [ml - (p + q) 2 ] 



(19) 



where Ai 



(m?— p 2 u) 



(p + q) 2 , u = 1 — u. In eqs.(18) and (19) we have rewritten the 



dispersion integral in terms of the variable u = (ml — p 2 )/(s — p 2 ) . 

Here we would like to make the following remark. As we noted earlier, the functions gi(u) 
and g2(u) represent twist r = 4 contributions to the two-particle photon wave function. To 
this accuracy, in eq.(19) we must take into account other twist r = 4 photon wave functions 
(see for example ||17|| ). Using the equation of motion, one can relate them to the three- 
particle wave functions of twist r = 4 with an additional gluon from heavy quark |T7|]. But, 
these three-particle wave function contributions, in general, are small and we will neglect 
them in further analysis. 

The remaining task is now to match eqs.(18) and (19) with the corresponding hadronic 
representation (see eq.(5)) and to extract the formfactors gi(p 2 ) and fi(p 2 )- As usual, in- 
voking duality, we assume that above certain threshold Sq = 35 GeV 2 (this value follows 
from two-point sum rules analysis) the spectral density p(s) associated with higher reso- 
nances and continuum states coincides with the spectral density from perturbative part. 
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This procedure is equivalent to writing (m 2 — P 2 )/(sq — p 2 ) in the lower limit of the integra- 
tion over u in eqs.(18) and (19) (for more detail see 0, jlSj). Finally applying the Borel 
transformation on the variable — (p + q) 2 — > M 2 to suppress both higher state resonances 
and higher Fock states in the full photon wave functions, we get the following sum rules for 
the formfactors: 



91 (P 2 ) 



+ 



+ 



-— e^ 7 

JB 



[ N c 
\l6vr 2 



1 du 



8 U z 



(ml — p 2 u 



1 - 



miu 



mi 



1 + 



mi 



p 2 u 



2e h 



2 2-^ 



+ 



m b (qq)e g I — 



u 



X(f>(u) - Am 2 b (9i ~ 92) 



u 2 M A 



LCdu 9 -^ 



A J 5 



u 



1 + 



p 2 + ml 
uM 2 



uM 2 



e b m b 



p 2 u 



(m 2 - p 2 ) 
2 «W 92 

2 9 
raj; — p z 



( m ft-p 2 ") 



+ 



r / 2\ ""6 



ra b m a 
— eM- 

JB 



\ 16tt 2 



5 rag — p 



1 - 



ragw 



ml — p 2 it^ 



x 



x ((. + et) =L^5 f ^? - V ) + (e, - e t ) m | f + / 



u 



11 



+ Amlp 2 ln 



miu 



+ m b e q (qq) 



1 dw 



5 u 



Am 2 , 

~ X4>{n) + ^(91- 92) + 



2 (p 2 + 
m 2 M 4 



-#2 



4 75 It 



+ 



-i + p 2 :™^ 



uM 2 



+ e b m b 



ml — p 2 



17 ? 7 



e 



(20) 



At the end of this section we give the result for the differential decay widths: 



dT 

ds 




v tb v; q 



m 5 B s(l 



\ 



ra 



m B s 



x 



\A\ a +\Bf 



H — |Cio | 2 



ra 



B 



f 2 (p 2 )+9 2 (p 2 ) 



(21) 
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where 



2 / 2 

s = p /m B , 

A = C e 9 ff g(p 2 )-2C 7 ^ gi ( P 2 ) , and 
B = C^f(p 2 )-2C 7 ^f 1 (p 2 ). 

4 Numerical Analysis 

For calculating formfactors fi{p 2 ) and gi{p 2 ) we use the following input parameters: 
m b = 4.7 GeV, s ~ 35 Ge^ 2 , / B = 140 MeF |§ 0, = 6«(1 - u) 0, 0. To the 
leading twist accuracy we use for g±(u) = 1 (see first reference in |TT[ ) and for gi(u) and 
#2 (w) the following expressions [p~3|] : 

0l(u) = -i(l-«)(3-«) (22) 
g 2 (u) = -\{l-u) 2 (23) 



The magnetic susceptibility \ was determined in ||I8|, x = —3.4 GeV 2 at the scale 



Hb ~ — ^6 , (qq) — —(0.26 GeV ) 3 . The Borel parameter M 2 has been varied in 

the region 8 GeV 2 < M 2 < 20 GeV 2 . Numerical analysis shows that the variation of M 2 in 
this region, changes the results by less than 8% . The predictions of the sum rules are very 
stable in this region of the Borel parameter and vary only a few percent with the changes 
of mft, s and fg within the intervals allowed by the two point sum rules for f B . 

The sum rules is reliable in the region m 2 — p 2 ~ a few GeV 2 , which is smaller than 
p 2 = ml- In order to extent our results to the whole region of p 2 we use some extrapolation 
formulas. We found that the best agreement is achieved by the dipole type formulas 



, 2N 3.74 GeV 2 , , 

9i(p 2 ) = _ p2 , (24) 

v rn f ' 

„ . 9x 0.68 GeV 2 , . 

MP 2 ) = _ p2 , (25) 

where m\ = 40.5 GeV 2 and m\ = 30 GeV 2 . For calculating differential and total decay 

widths, we need the values of Cg^ , CV and C\q coefficients and the explicit forms of the 
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formf actors g(p 2 ) and f{p 2 )- These formfactors are calculated in \TT\: 

1 GeV 



9(P 



f(P 2 



'1 _ P \2 
5.6 2 ^ 



0.8 GeV 

y L 6.52 1 



(26) 
(27) 



The values of the Wilson coefficients C7 and C±o are taken from & 1C] as 



CV = -0.315 , C 



10 



-4.642 , 



and the expression Cg^ for b —>■ s transition, in the next-to-leading order approximation is 



given as (see [|T!|) 



r< e ff 



C 9 + 0.124w(s) + g(rfi c , s)(3d + C 2 + 3C 3 + C 4 + 3C 5 + C 6 
~g(m q , s)(C 3 + 3C 4 ) - ^g(m b , s)(AC 3 + AC 4 + 3C 5 + C 6 ) - 
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(3C 3 + C 4 + 3C 5 + C 6 ) , 



(2f 



with 



-0.249 , C 2 = 1.108 , C 3 = 1.112 x 10" 2 , C 4 = -2.569 x 10" 



C h = 7.4 x 10 



-3 



Ca 



-3.144 x 10" 2 , Cg = 4.227. 

The value of Cg for b d transition, can be obtained by adding to eq.(28) the term 
X u [g{m c , s) - g(m d , §)) (3 C x + C 2 ), where 

v ubV* d 



For obtaining these values we used Aq C d = 225MeV, sin 2 9w = 0.23, m 4 = 176 Ge^, raw = 
80.2 GeV and m q = m^/m^. In the above formula w(s) represents the one-gluon correction 



to the matrix element Og and explicit expression can be found in |]10| , while the function 
g(rfi q , s) arises from the one loop contributions of the four quark operators 0\ - Oq (see for 
example || [K|), i-e. 



g(m q ,s') = --lnm q + — + -y q --(2 + y q )^Jll-y q + +^Q(l-y q )x 

x (/n ) + V i^ - m ) + 6(1/, - l)arctg-= L— ) (29) 
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with y q = m q 2 /s' } and s' = y 2 jm\. 

For a more complete analysis of the above decay, one has to take into account the long 
distance contributions. In the case of the J/ip family, this is accomplished by introducing 
a Breit-Wigner formula through the replacement (see |20| ) 

3tt myBrjV -> l+l')?} ' 

where rdy = vnyjra^^ f tot = T/rrib . The masses and decay widths of the corresponding 



g(rh c ,s')^g(rh c ,s')-- £ / \ o . . (30) 

a v=J/M' s'-m v +im v Tl 



mesons are listed in |^TJ. In Fig. 3 we present the differential decay rate for B s — ► yU + /i~7 
decay (behavior of the differential decay rate for B s — > e + e~7 decay is similar) as a function 
of s, with and without resonance (J/ip and if)') contributions. From this figure we see that 
the contribution from soft photons, corresponding to large s region is negligible. 

Using the above mentioned values of the parameters and | V^t, V^.* | = 0.045 , |Vtf,V^| = 0.01, 



r(B s ) = 1.34 x 10 12 s , r(B d ) = 1.5 x 10 12 s [pTfl , for branching ratios we get (without 
the long distance contributions): 



B(B S -> e + e"7) = 2.35 x 10~ 9 

B(B S -> /i + /i~7) = 1.9 x 10~ 9 

B(B d -> e + e~j) = 1.5 x 10~ 10 

B(B d -> // + /i~7) = 1.2 x 10~ 10 (31) 

For comparison we present also the constituent model prediction (at /# = 140 MeV , 
m s = 0.57 GeV , m d = 0.35 GeV) [§: 



B{B S - 


-> e + e 


7) 


= 3 x 10~ 9 


B(B S - 




7) 


= 2.3 x 10" 9 


B{B d - 


-> e + e" 


7) 


= 4 x 10~ 10 


B{B d - 




7) 


= 3 x 10~ 10 



(32) 

We see that the constituent quark model and light cone sum rules method predictions 
on the branching ratios are very close. Let us compare our results on branching ratios with 
those of pure leptonic decays. The rates for the pure leptonic decays are (see for example 
11) 

V" "' IT \ -:|-|2/-.2 

167T 3 



a 2 G%.f% m B mf , , 

r(5, - /+r) = 9 > ^*/Cf (33) 
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If we use the value of fs s — fs d — 140 MeV, for the corresponding Branching ratios we 
get: 



B{B S - 


-> e + e 


") 


= 3 x 1(T 14 


B{B S - 




') 


= 1.3 x 10"' 


B(B d - 


-> e + e~ 


") 


= 2.1 x 10~ 


B(B d - 







= 9 x 10 -11 



From these values and eq.(30) it follows that, the radiative decays dominate over the pure 
leptonic decays in the corresponding channels and B s — > e + e~7 decay mode has a larger 
branching ratio. Few words about the experimental detectabilty of these processes is in 
order. In future 5-factories and LHC approximately 6 x 10 11 (2 x 10 11 ) B d (B s ) mesons are 
expected per year. Therefore the decays B s ^ — > Z + Z~7 are expected to be quite detectable 
in these machines. 

In conclusion, we have analyzed the rare B q — > Z + Z~7 decays in SM and obtain the 
branching ratios for B s — > Z + Z~7 to be around 2 x 10~ 9 and around 2 x 10~ 10 for B d — > Z + Z~7. 
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Figure Captions 

1. Feynman diagrams in the Standard Model for bq — > 

2. Diagrams describing the perturbative and non-perturbative contributions to the corre- 
lator function (4). 

3. Differential decay rates of B s — > yU + /i~7 versus s = p 2 /m 2 B . Here the thick line corre- 
sponds to the case without the J/ip, ip' and the thin line with the J/ip, ip' contributions. 
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